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SUMMARY 


A calculation method is presented for sound propagation over an impedance 
discontinuity in flat ground with a homogeneous, still atmosphere. The method is 
based on an approximate solution to a two dimensional boundary integral equation 
formulation of the problem, which expresses the wave field as the solution for 
homogeneous ground plus an integral over half of the boundary. Through recognising 
this integral as a generalised Fourier integral, asymptotic methods are applied to 
evaluate the part of the integral most expensive to compute by numerical quadrature. 
Single frequency excess attenuation results for propagation from a point source 
above rigid ground to a receiver above absorbing ground are discussed. The results 
are applied, with air attenuation and A-weighting, to a notional jet engine noise 
source; simple trends are noted. 


INTRODUCTION 


The problem discussed in this paper is propagation from a point source in a 
homogeneous still atmosphere above flat locally reacting ground. Efficient 
calculation methods for the wave field above acoustically homogeneous ground are 
well known (e.g. ref. 1). More recently sound propagation over impedance 
inhomogeneities has been theoretically examined; a thorough review is given in 
reference 2. A limitation of the accurate calculation methods is their 
computational expense . 

Here we focus on propagation over a single straight line impedance discontinuity 
which lies perpendicular to the direct source-receiver propagation path. A 
development to an existing calculation method is described which significantly 
reduces the computational expense. 

The improved calculation method is applied to grazing incidence propagation from 
a source above a rigid surface to a distant receiver above absorbing ground. 
Monofrequency excess attenuation results are examined and some simple trends are 
observed. The results for a l-5m high receiver are applied, with air attenuation 
and A-weighting, to a notional jet engine noise source at l*5m height. Again some 
simple trends are noted. 
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CALCULATION METHOD 


t 


Description of the Problem 

Figure 1 illustrates the problem. A point source with harmonic time dependence 
( e ~io)ty i s situated over a flat locally reacting surface of infinite extent. The 
surface is divided by a straight line into two half planes. Each half plane is 
acoustically homogeneous and characterised by a frequency dependent complex 
admittance (the inverse of the normalised acoustic surface impedance). We are 
interested in evaluating the acoustic potential at a point in a vertical half plane 
that is bounded by the surface, passes through the source and is perpendicular to 
the line of the admittance discontinuity. For the mathematical description we will 
use right-handed Cartesian coordinates Oxyz as indicated in Figure 1, the y-axis 
vertical and the surface in the plane y=0 . The source and receiver coordinates are 
(0,h s ,0) and (L,h r ,0) respectively. The admittance discontinuity is along the line 
x=X in the surface . 


An Existing Calculation Method 


First we consider the related problem in which the source is replaced by an 
infinitely long coherent line source, parallel to the admittance discontinuity. 
This cylindrical wave propagation problem is mathematically equivalent to the two 
dimensional problem which is illustrated in Figure 2. From the mathematical 
expression of this problem as a two dimensional boundary value problem (consisting 
of the Helmholtz equation and suitable boundary conditions) the following boundary 
integral equation can be derived (ref. 3): 


- G ( 32 (t 1 ,t 2 ) + ik((5,-p 2 )J* oo *(s, t, )G 0 2 (s, t 2 )dx. 


( 1 ) 


In this equation ^=*(0,1^) is the source position, t 2 =(L,h r ) is the receiver 
position, and /3 1 and (3 2 are the admittances of the two halves of the boundary. The 
integration is over the interval ^=(-00 ,X] ; this is the part of the boundary with 
admittance /3 1 . s=(x,0) is a point in the boundary. For two points a and b, <t>(a,b) 
denotes the acoustic potential detected by a receiver at a when insonified by a unit 
source at b; G^(a,b) (where |3= = /3 1 or (3 2 ) denotes the same quantity in the simple 
case when the boundary has homogeneous admittance j3. Efficient methods for 
evaluating the solution in this simpler case have already been developed (refs. 4,5). 


Equation (1), which describes an inhomogeneous admittance boundary problem, can be 
solved accurately for by the boundary element method (refs. 3,4,6). We 

consider here an approximate but less computationally expensive method of solution. 
To develop this we make the physically plausible assumption that the potential in y 
is what it would be if the whole boundary had admittance |3 1 (refs. 7,8). Thus 
$(§,£i) equation (1) is replaced by G^^s,^), giving the following approximation 

to <fc(t 2 , t, ) : 


where 


4 > a(£ 2 ,£i) = G 02 <£i ,t 2 ) 


+ i(0,-0 2 )I(X) 


•x 


I(X)= k 


G /S i < S . t n 2 ( § , t 2 ) dx • 

—00 


( 2 ) 


This approximation avoids using the boundary element method. 
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Using < 1 , a(£2>£i) we can calculate Qa(£ 2 >£,)> an approximate cylindrical wave 
reflection coefficient. Let d and D denote the distances from source and from image 
source (at (0,-h s )) to receiver, respectively. Assuming a source with unit volume 
flow rate amplitude, 

*a(£ 2 ,£i> - -^HjO(kd) + Q A (£ 2 ,£ 1 )[-jH^)(k D )] , 

where H^ 1 ) is the Hankel function of the first kind of order zero. 

For propagation over short distances, the difference 201og 10 itf’A^ ,t )l 
-201og^Ql < * ) (£ 2 .£ 1 ) I . where ^(tj.t^ is calculated by the boundary element 1 method, has 
been found to be around 0-ldB (ref. 6). This suggests that is an accurate 

approximation to Q(t 2 ,t,), the exact cylindrical wave reflection coefficient for the 
two dimensional problem illustrated in Figure 2. Also, it has been argued that 
(ref. 6), for a receiver in the far field of the image source (kD>l) , Q(t 2 ,t ) is an 
accurate approximation to the spherical wave reflection coefficent, q, for the three 
dimensional problem illustrated in Figure 1. If the point source in the three 
dimensional problem has unit volume flow rate amplitude, then the acoustic potential 
at the receiver position is 

e ikd e ikD 
^ 4ird ^ 4irD 

Replacing q with Q(t 2 ,t,), which is approximated by Qa(£ 2 .£,)> we obtain an 
approximation for <p : 

e ikd e ikD 

*A = - -2£d " Q a(£2,£,)-4^d • (3) 

We assume throughout the rest of the paper, without further comment, that this is a 
good approximation for <p. 

The main computational expense in this approximate calculation method is in 
evaluating the integral I(X). In previous calculations (ref. 6) I(X) was evaluated 
numerically after first replacing the lower limit of integration, - 00 , by a 

^• c ^- an ^'3.y large negative value. Unfortunately the integrand in equation (2) is 
usually highly oscillatory over the range of integration, making numerical 
integration an expensive process. Here we derive a semi-analytical method of 
evaluation which deals efficiently with the part of the integral that is most 
expensive to evaluate numerically. 


The Improved Calculation Method 

We begin by examining the general behaviour of the integrand in equation (2). 

We note that, for a receiver at the point s in the boundary, we can write 

4iG /S(§,t) - 1 ) (k| t-s | )R|j(s , t) , 

where Rp(s , t)=l+Q^(s , t) , and Q^(s,t) is the cylindrical wave reflection coefficient 
for a homogeneous surface of admittance 0. As x is increased from - 00 , the real and 
imaginary parts of H^O(k|t-s|) oscillate in a well defined fashion, while Ro(s,t) 
changes less rapidly. In fact, from the asymptotic expansion of the Hankel function 
at large argument, we see that if we factorise, 

G /3 ( s,t) = eik|t-si s(§i t) , 
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then S(s,t) is smooth as a function of x compared to ei k| £ - § ( when k|t-s| is large. 
This observation suggests that I(X) can be usefully written in the form of a 
generalised Fourier integral: 


I(X) =[ X f(x)e ik S( x )dx , 
** —00 


( 4 ) 


f(x) = kGp, (s,t,)G|5 2 (s,t 2 )e i k g( x ) , 
g (x) = g, (x) + g 2 (x) , 

g,(x) = |t,-S|, g 2 (x) = lt 2 -s|. 

Notice that g(x) is the distance from source to receiver via the point s. The 
location on the boundary of the geometrical reflection point, x r , is therefore given 
by g'(x r )= 0 . When s is sufficiently distant from t 1 , t 2 , and (x r , 0 ), f(x) is a 
slowly changing function of x compared to eikg(x) . This fact allows us to use 
simple asymptotic methods to help evaluate I(X). 


where 


and 


To introduce the asymptotic analysis we consider first what proves to be the 
simplest type of configuration to deal with. This has x r well outside 7, and and 
t 2 at least one wavelength from the boundary. For this type of configuration we may 
integrate I(X) by parts to give 

I(X) = J,(X) + R,(X) ( 5 ) 

where 


J,(X) = fm e lk g( x ) 


ikg'(X) 


R,(X) 


[ x d_ 

f(x) ^ 

dx 

—00 

ikg' (x)J 


: ikg(x) dx 


and then integrate R, (X) by parts to give 

R,(X)-J 2 (X)+R 2 (X) . 

where 


J 2 (X) = 


g"(X) f'(x)l f(X)e ik e( x ) 

g’lxy f (X) J ( ikg’ (X) ) 2 


R 2 (X) 


f X d_ 

d 

f (x) 1 

1 

dx 

—00 

dx 

ikg’ (x) 

ikg' (x) 


e ikg(x) dx 


If f(x) were completely independent of k, we could apply the Riemann Lebesgue 
lemma to show that, for n~l , 2 , 

R n (X) = o(k~ n ) , k -> 00 . 

In fact f(x) depends weakly on k, but f(x) approaches a limit independent of k as 
k->oo with other variables fixed (ref. 2 , p. 585 ). Thus I(X) has the following 
asymptotic approximation: 

I (X) - J , (X) + J 2 (X) , k ^ 00 . 
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J 1 (X) and J 2 (X) are the first and second terms in what is approximately an 
asymptotic expansion of I(X) in inverse powers of k. When k is large enough, 
J 1 (X)^J 2 (X), so that we can safely approximate 

I (X) * J,(X). (6) 

The above arguments do not tell us how large k should be in any particular case 
for approximation (6) to be valid. However, it is plausible that R 2 (X)<J 2 (X) when 
^ 2 (X) i (X) . Thus we can estimate the relative error in approximation (6) by the 
following upper bound on i J 2 (X)/J 1 (X) | : 

E r (X) 


f g" (X) 3 

r_L_ + 1 1 

1 1 

llg’ (X) I 21 

Lg, (X) g 2 (X)J 

Jklg' (X) 1 


(To obtain this expression, |f*(X)/f(X)| has been replaced by 

(3/2) [l/g 1 (X) + l/g 2 (X)], which is expected to be an upper bound on |f f (X)/f(X)| in 
all cases (ref. 2, p.598).) We can estimate the absolute error in approximation (6) 
by the following upper bound on | J 2 (X) | : 

lf(X)| 


E a (X) - E r (X) 


klg f (X) | 


Both E r (X) and E a (X) are infinite at X=x r and tend to zero as X-»-oo. Moreover, a 
graphical examination of E r (x) and E a (x) suggests that they are monotonic in 
(-°o,x r ) , for typical geometries, admittance values and frequencies. 


We move on to consider configurations for which still X<x r but x r -X is small 
enough for E r (X) and/or E a (X) to be unacceptably large. For the moment we require 
that both the source and receiver are many wavelengths above the boundary. The 
following breakdown of the integral is used: 

I (X) - I (t ) + K , (7) 


where K is the integral over a truncated interval 7T=U>X], 

f X 

K - kj ^/? 1 (§>J 1 )G|3 2 (s,t 2 )dx , 


which will be evaluated numerically. We will choose r so that we can satisfactorily 
approximate 

KO - J,(t) . (8) 


To reduce the expense in evaluating K numerically we want to choose r as close 
to x r as possible while still insisting that approximation (8) should satisfy 
certain relative and absolute error criteria. We can uniquely define two upper 
limits, t r and r a , for r by 

E r (r O — 0 * 2 , 


E a (r a ) - e , (9) 

where e is an arbitrary positive constant. r<r r ensures that J 2 (r) is sufficiently 
small compared to J n (0 for E a (r) to be an accurate estimate of the absolute error 
in approximation (8). Therefore, if also 7 < 7 a , then the absolute error in the 
approximation (8) is <e . Thus by taking r to be the minimum of r r and r a we ensure 
that the error made in replacing I(r) by J^r) in equation (7) is <c . 
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We can now consider the more general configuration in which x r may be anywhere 

in relation to X, but the source and receiver remain many wavelengths above the 

boundary. If X is less than x r , one of the above calculation methods applies. If X 
is greater than x r , reciprocity can be invoked (reflect the problem in the plane 
x=L/2, then swap the source and receiver) and then one of the above methods applied. 

Finally we note why the source and the receiver have so far been kept at least 

one wavelength above the boundary. If t 1 (or t 2 ) is very close to y the 

approximation (6), which involves neglecting the integral R 1 (X) in equation (5), 
breaks down. This is because, for x in a small range of y around t 1 (or t 2 ) , f(x) 
changes rapidly with x. Thus the derivative of f(x) in the integrand of R, (X) is 
very large . 

To avoid the consequent inaccuracies which may occur when t n or t 2 is within one 
wavelength of y, an additional criterion is used for the choice of r . We require r 
to be small enough so that the line x<r in y is always at least one wavelength from 
1 1 and 1 2 . 


GRAZING INCIDENCE RESULTS 


Monofrequency Excess Attenuation 

We can use the method described above for estimating ^ to examine propagation 
over flat ground through a homogeneous still atmosphere. The monofrequency excess 
attenuation over geometrical spreading due to the presence of the ground can be 
approximated by 

ikd 

Al - 201og, 0 | 5 S a -^-| d» . (10) 

Propagation from a source above rigid ground to a receiver above absorbing 
ground has been examined. To model this problem |3 1 was set to zero and the 
dependence of (3 2 on frequency was calculated by the Delany and Bazley semi-empirical 
formula (refs. 9,10), with an effective flow resistivity of 10^kgs -1 m” 3 . This value 
was chosen as being a low value for grassland (ref. 11) . (It is found that using 
two or three times this flow resistivity value causes only a small reduction in the 
magnitude of the A^ results, and no change in the trends was observed.) Six 
configurations of source and receiver heights (h s and h r ) were examined. 

Specifically, heights of 5m, l*5m, and 0-5m were used, with h s >h r in all cases. For 
each h s and h r combination calculations were carried out at four distances: L=250m, 
500m, 1km, and 2km. 

We will examine the significance of the proportion of rigid ground between the 
source and receiver. We can define a useful variable, p r , by 

p r = X/L . 

When the impedance discontinuity is between the source and receiver (0<X<L) , p r 
gives the proportion of rigid ground between the source and receiver. 

A sample of the results examined is shown in Figures 3a and 3b. Plots like 
those shown were calculated for all the octave band centre frequencies between 100Hz 
and 5kHz. Notice that, as is of course expected intuitively, when p^XO or p r >l the 
modelled ground behaves as an acoustically homogeneous plan, absorbing or rigid, 
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respectively. We therefore now concentrate the investigation only on the ranee 
0<p r <l. 6 

Figure 4(a) illustrates a simple curve shape that occurs whenever both kh s and 
kh r are small enough. Half of the plots examined were of this type. We see that Ai 
increases monotonically with p r . The gradient of each curve is greatest at the ends 
of the range of interest, i.e. near p r =0 and p r =l. At low enough frequencies the 
curves straighten out. 

Figure 4(b) illustrates the disintegration of the orderly patterns seen in 
Figure 4(a) that occurs when h s , h r , or the frequency is increased sufficiently. 
About a third of the plots examined showed this type of disorderly pattern. 

Figure 4(c) illustrates a different pattern that sometimes occurs when the 
source is higher than the receiver, but neither are so high above the boundary that 
the disorderly pattern seen in Figure 4(b) occurs. In Figure 4(c), A 1 is less 
dependent on the location of the impedance discontinuity when p r <0-3. Notice that 
the right half of this plot shows the features observed in Figure 4(a). About a 
sixth of the plots examined showed this pattern. 

When the orderly patterns seen in Figure 4(a) and the right hand side of Figure 
4(c) occur, there is usually a range of octave band centre frequencies around 1kHz 
at which some or all of the curves on a plot are separated, in most of the range 
0<p r <l, by about 3dB. This approximate 3dB increase in A^ per doubling of L when 
propagation is over an admittance discontinuity occurs only for a range of values of 
kjL- The start of this inhomogeneous absorbing ground effect corresponds with the 
start of a 6dB separation of the A x curves at p r <0, which occurs when kL=2300. This 
6dB increase in A]_ per doubling of L is a homogeneous absorbing ground effect which 
has been predicted theoretically (ref. 12). Figure 5 illustrates these 
observations . The range of plots examined show that the inhomogeneous ground effect 
fails to occur when kh s or kh r is large. This failure is observed in the left hand 
side of Figure 4(c) . 


Jet Engine Noise 

We move on to examine the excess attenuation of a notional broad band 
environmental noise source. A simple spectral shape representative of a jet engine 
Full thrust is chosen. The free field lm third octave band sound pressure level 
is taken as constant up to 200Hz, above which frequency it is reduced by 0-8dB per 
third octave band. Third octave band excess attenuations due to the presence of the 
ground are approximated here by A]_ values given by equation (10), using the band 
centre frequencies. To make the calculation more realistic we include the B.S.5727 
(1979) third octave band free field air attenuations for 20 *C and 70% relative 
humidity . 

We consider only source and receiver heights of l-5m so that the simple A^ 
pattern illustrated in Figure 4(a) dominates the results. The process of intensity 
addition over the third octave bands will produce more moderate excess attenuations 
for the broad band noise than those calculated for monofrequency sound. 

The excess attenuation of the broad band source noise caused by the presence of 
the ground and by air absorption along the propagation path is 

A 2 “ Sx - (S 2 + 201og, 0 L) dB(A) , 
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where is the total A-weighted sound pressure level of the notional source at lm 
distance in the free field, and S 2 is the total A-weighted sound pressure level that 
we calculate at the receiver position. 

Figure 6 shows how A 2 depends on p r and L. We can see that it is not possible 
to predict the A 2 values at intermediate p r from a linear interpolation between the 
values at p r =0 and p r =l . Notice however that, in the range 0 * 25<p r <0 * 75 , A 2 appears 
to vary linearly with p r . Also, in this range, A 2 increases by 5dB(A) per doubling 
of L. This dependence on L is in reasonable agreement with a commonly used 4dB(A) 
extra attenuation (of perceived noise level) per doubling of receiver distance from 
an aircraft source very close to the ground (of unspecified admittance) (ref. 13). 

At L=lkm and 2km, A 2 varies linearly with p r in the range 0<p r <0-75. This is 
useful because it means that A 2 (p r ) can be estimated from A 2 (0), the value for 
homogeneous absorbing ground, which is easier to calculate. The simple predictive 
equation, which is shown in Figure 5 for L=lkm and 2km, is 

A 2 (p r ) “ A 2 (0) - 14 • 7p r dB(A) , 0 < p r < 0-75, (11) 

Unfortunately we know of no practical results with which to compare this equation. 


CONCLUSIONS 


An improved calculation method has been presented for sound propagation over a 
straight line impedance discontinuity in flat ground. The method is restricted to 
the case when the impedance discontinuity is perpendicular to the direct source to 
receiver propagation path. The method is derived from an asymptotic analysis at 
large wavenumber of an approximate solution of a two dimensional boundary integral 
equation. Accuracy is adequate for the purpose of examining environmental noise 
propagation in ideal conditions. A limitation of the method is the assumption of 
homogeneous still air and flat ground. 

Results for long distance grazing incidence monofrequency propagation show that 
the dimensionless heights (height multiplied by wavenumber) of the source and 
receiver above the ground are as important as the location of the impedance 
discontinuity. When these dimensionless heights are small enough, the results are 
very orderly, as illustrated by Figures 4(a) and (c) . 

Theoretical results for the excess attenuation, including air absorption, of a 
broad-band A-weighted notional environmental noise source have been examined. A few 
simple trends have been noted, in particular, equation (11). 
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Fig.l The problem in three dimensions. 
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Fig. 4a Monofrequency excess attenuation over geometrical spreading at four long 
distances, plotted against the proportion of rigid ground. h s “h r «l*5m, f“250Hz. 
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Fig. 4b Monofrequency excess attenuation over geometrical spreading at four long 
distances, plotted against the proportion of rigid ground. h s =5m, h =1.5m, f=2kHz 













